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Abstract. - We study a 3-dimensional Ising model in which the tendency to order due to 
short-range ferromagnetic interactions is frustrated by competing long-range (Coulombic) inter- 
actions. Complete ferromagnetic ordering is impossible for any nonzero value of the frustration 
parameter, but the system displays a variety of phases characterized by periodically modulated 
structures. We have performed extensive Monte-Carlo simulations which provide strong evi- 
dence that the microphase separation transition between paramagnetic and modulated phases 
is a fluctuation-induced first-order transition. Additional transitions to various commensurate 
phases may also occur when further lowering the temperature. 



In a variety of systems a tendency to order induced by short-range interactions is frus- 
trated by long-range competing effects. In diblock copolymer melts, cross-linked polymer 
mixtures and interpenetrating networks, short-range forces between mutually incompatible 
components would drive a phase transition at low enough temperature, but total segrega- 
tion is forbidden by constraints due to covalent bonds or crosslinks. Instead, a microphase 
separation transition may occur, where the system forms phases with periodically modulated 
structures such as lamellar, hexagonal, or cubic phases^, |^, |^j. Similarly, self-assembly in 
water-oil-surfactant mixtures results from the competition between the tendency of oil and 
water to phase separate and the stoichiometric constraints due to the presence of the surfactant 
molecules that act as the electroneutrality condition in a system of charged particles^, |j. In 
supercooled liquids, the dramatic slowing down of the structural relaxation which leads to 
the glass transition has been ascribed to the presence of frustration-limited domains whose 
formation comes from the inability of the locally preferred arrangement of the molecules in 
the liquid to tile space periodically^. Others examples include frustrated phase transition 
in doped antiferromagnets[Q and metal-ammonia solutions^), or else pattern formation in 
various nonequilibrium systems^. 

The physics of the above mentioned situations appears to be reasonably well described by 
lattice or continuum models in which both short-range ordering and long-range Coulombic 
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frustrating interactions are present. In this letter, we consider a Coulomb frustrated Ising 
fcrromagnet whose Hamiltonian is given by 

H=-jj2s iSj + ^j:^i, (i) 

<i,j> ijtj %3 

where J 7 Q > 0, Si = ±1, < i,j > indicates a sum restricted to nearest neighbors, and r»j 
is the distance between the sites i and j on a three-dimensional cubic lattice. The two main 
questions that we address here by combining analytical results and Monte-Carlo simulations 
are the following: (i) what is the nature of the microphase separation transition from the 
paramagnetic to the "modulated" phases and how does it vary with varying frustration, i.e. 
varying ratio Q/J? (ii) What is the nature of the modulated phases? 

The mean- field approximation to the Coulomb frustrated model is easily derived ||. We 
introduce t^{Q) as the Fourier transform of the full pair interaction potential in Eq.Q e.g., for 
small wave vectors and a 3 — d cubic lattice, Tk(Q) is approximatively given by 

r k (Q) ~ -6 + A(k) + |g (l + 6A(k) + C ^^M + ..^ (2 ) 

where A(k) = Yl a =x,y z = 2 J2 a =x, y . z ( l ~ cos(k a )) is the Fourier transform of the 

lattice Laplacian, b and c are constants, and J has been set equal to 1. For Q ^ 0, 7k(Q) 
attains its minimum, r m (Q) = Mm^{T^(Q)}, for a set of nonzero values of the wave vector, 
{k TO (Q)}. The mean-field treatment predicts the occurrence with decreasing temperature of a 
second-order phase transition at T C (Q) = —T m (Q); T C (Q) goes continuously to T®, the critical 
temperature of the unfrustrated ferromagnet, as T% — T C (Q) ~ Q 1 / 2 when Q — ► 0. At the 
transition, the ordering is characterized by a nonzero wave vector, corresponding to one of the 
k m (Q)'s, which for small enough values of Q is zero in two directions and behaves as Q 1 ^ 4 in 
the third one. The mean spherical version of Eq.(Q), in which the spins are taken to be real 
numbers with the global constraint that their mean square value is equal to one, leads to a 
quite different behavior^ The most significant feature is the existence of an "avoided" 
critical point at T® |Tc| ] . The microphase separation transition is again of second-order and the 
ordering wave vector is one of the {k TO (Q)}'s; but the critical temperature is now given by 
T~ 1 (Q) = J d 3 fc/(27r) 3 (Yk(Q) — T m (Q)), and when Q — > 0+ it goes to a temperature which 
is much less than T c °. Actually, within the mean spherical treatment, the transition is very 
sensitive to the form of the long range interaction: whereas a 1/r interaction leads to T C (Q) > 0, 
a slightly different form obtained by retaining only the leading 1/A(k) term in Eq.(j^) drives 
T C (Q) to zero for Q < 16|jl^|. This latter feature results from the fact that the minimizing wave 
vectors k m (Q) form a 2 — d surface in k— space (they form a finite set for the 1/r interaction). 
A similar feature is encountered in the Landau-Ginzburg- Wilson Hamiltonian that describes 
symmetric diblock copolymers . It has been argued by Brazovskii[^T|, on the basis of a self- 
consistent Hartree approximation, that such systems display a fluctuation-induced first-order 
transition that takes place at a nonzero temperature Tq{Q), with T° — Tq(Q) ~ Q 1 ^ 4 when 
Q — > 0+. Additional support to this prediction has been given by renormalization group 
analyzes |]l2f . 

In order to check the relevance of these various predictions, we have performed extensive 
Monte-Carlo simulations of the 3— dimensional Coulomb frustrated Ising ferromagnet. Our re- 
sults provide, we believe for the first time, direct numerical evidence supporting the Brazovskii 
scenario of a fluctuation-induced first-order transition for such a model. The simulations have 
been done for many values of the frustration parameter Q. Long-range interactions lead to 
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very time-consuming runs and we have considered lattices of L? sites for L ranging from 4 to 
24 with periodic boundary conditions (anisotropic lattices of various sizes up to 12 x 12 x 30 
have been also used). To treat properly the Coulomb interactions Ewald sums have been 
used; the site-site pair terms are calculated once for all at the beginning of the run and are 
stored in an array for the entire run: therefore, a large number of reciprocal vectors can 
be included in the Ewald sum to ensure a very good accuracy for the calculation of the 
Coulomb potential. The constraint of zero total magnetization has been enforced by always 
considering the simultaneous flip of a pair of up and down spins (note that the presence of 
long-range interactions and the relatively fast equilibration for system sizes studied reduce the 
advantage of considering cluster algorithms). Each run consists then of 2 — 6 x 10 4 Monte-Carlo 
Sweeps (MCS) per lattice site for equilibration followed by 6 x 10 4 — 10 5 MCS to calculate the 
equilibrium quantities. In order to verify that the equilibrium state is indeed reached and that 
the results do not depend upon the initial configuration, each run is performed by starting 
both from a high-temperature configuration and from the ground state. 

In Fig. 1, we display the microphase separation transition from paramagnetic to modulated 
phases in the (scaled) temperature-frustration diagram for a dozen of values of Q. Also shown 
on the same diagram are the mean-field and mean-spherical predictions. Transitions associated 
with very small values of Q are of course inaccessible to simulations in finite-size boxes (at 
least two lamellae characterized by opposite magnetizations must be present in the simulation 
box to satisfy the constraint of zero total magnetization, which limits the width of the lamellae 
and by consequence the smallness of the frustration); however, the results are compatible with 
a power-law behavior, 

T c ° - T (Q) ~ Q K (3) 

with k ~ 0.25 — 0.35, and the period of the lamellar structures is in good agreement with 
27r/|k m (<2)|. As stressed before, the transition is first-order. To locate the transition and 
determine its order, we have monitored as a function of both decreasing and increasing 
temperature the energy histograms, which provides the total energy per spin E of the system, 
the heat capacity C v =< (E— < E >) 2 > and the Binder cumulants; we have also studied 
the (potential) order parameter \M^\ = | < > |, where is the Fourier transform of the 
spin variable along the three directions of the lattice. Results for Q — 0.4 are shown in Fig. 
2. A hysteresis loop between heating and cooling runs is observed in E versus T (Fig. 2a), 
and the distribution of the energy histograms around the transition temperature To exhibits 
a forbidden region for the energy (Fig. 2b); they both indicate a first-order transition. Also 
supporting this conclusion is the apparent jump in the order parameter and the strong system 
size-dependence of the heat capacity (e.g. MaxTCt,(T) ~ LP with p > 1). Similar behavior is 
observed for the other values of Q. However, the transition becomes more and more weakly 
first-order as the frustration decreases. 

To provide a more complete investigation of the phase diagram, we have also studied 
analytically and numerically the T — properties of the system. In doing so, we have 
assumed, as indeed observed in our Monte-Carlo simulations and also in several related two- 
dimensional models Jl3|, |l4| and in the three-dimensional Ising model with dipolar competing 
interactions [l5|], that the ground states consist in periodically modulated structures. When 
Q — > +oo the ground state is an anti-ferromagnetic (Neel) state (lxlxl) and it is stable 
for Q > 15.3. For 15.3 > Q > 14.6, the system passes though a succession of "anisotropic 
checkerboard" phases of the form (1 x 1 x n) with n a positive integer. A "tubular" phase 
(1 x 1 x oo ) is stable for 14.6 > Q > 6.2, and in a narrow range of Q (6.2 > Q > 5.2), one obtains 
a second infinite sequence of (1 x n x oo) phases. For smaller values of Q, the ground state 
is found by lamellar phases whose width increases with decreasing Q. We have checked that 
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Fig. 1. - Location of the microphase separation transition in the T/T° — ln(Q) diagram. (Recall 
that J = 1.) The dots denote the Monte-Carlo results (T° ~ 4.51) The full line corresponds to the 
mean-field prediction (T° = 6) and the dotted line to the mean-spherical prediction (T° = 3.91...). 



more complicated phases involving mixtures of lamellae with different widths are not stable 
at zero temperature. When Q — > and for the purpose of calculating the energy, the Fourier 
transform of the Coulomb potential can be well approximated by 1/A(k). For a lamellar 
phase of period 2m in the z— direction, the only nonzero values of |Sk| are for k = (0, 0k z ) 
and \Sk z \ = l/(msin(fc z /2)) with k z = 7r(2n + l)/m and n = 0, ...,m — 1. Therefore, after 
some calculation, one obtains an explicit expression for the corresponding energy per spin, 

E(m) = - J(3 - -) + Q (lm 2 + J) , (4) 
m V o A / 

showing that the lamellar width m goes as (nQ/6) 1 ^ 3 when Q — > 0. This dependence is 
analogous to that derived for symmetric diblock copolymer systems in the so-called strong 
segregation limit 0, w], whereas the behavior at the transition line to the paramagnetic state 
rather corresponds to the weak-segregation limit result [jj], || . 

At nonzero temperatures, one may expect that the system undergoes a series of transitions 
to various commensurate, and at temperatures close to T a (Q) incommensurate, phases just 
like in the ANNNI model [^6|. Such a behavior is indeed observed in the mean- field theory §, 
but a sequence of modulated phases is hard to detect in Monte-Carlo simulations of finite- 
size systems. We have focused our study on Q = 0.144 for which the ground state is a 
lamellar structure of half-period 2 and the ordering at the microphase separation transition is 
characterized by a modulation half-period close to 2.5. (Other modulations could also occur 
at intermediate temperatures.) To observe at least these two types of order in a simulation 
where only restricted values of the wave vectors are allowed, we have used anisotropic lattices. 
For instance, with a 12 x 12 x 24 lattice, phases with modulation of half-period of 2.4 are 
allowed in the z-direction. Figure 3 illustrates at the level of the order parameter |Mk| the 
passage from the low-temperature state formed by periodic lamellae of width 2 (k z — tt/2) to 
a mixed or incommensurate phase characterized by an average half-period of 2.4 {k z = 5ir/12) 
for 1.85 > T > 1.75. For T > 1.85, the system is disordered and all components of \M^\ are 
zero. Corresponding spin configurations are displayed in Fig. 4. 

In this letter, we have studied the properties of the microphase separation transition and 
of the modulated phases observed in a three dimensional Coulomb frustrated Ising model. 
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Fig. 2. - Monte-Carlo results for (a) the total energy per spin E as a function of T (open and filled 
circles correspond to heating and cooling runs, respectively); (b) energy histograms around T — 1.5 
(full line in (a)) for a frustration parameter Q = 0.4. 
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Fig. 3. - Variation with T of the |M k | with k = (0,0, A:*) for (a) k z = tt/2 and (b) k z = 5tt/12 
from the Monte-Carlo simulations (for Q = 0.144 and a box size =12 x 12 x 24). The microphase 
separation transition occurs for To ~ 1.84 and it is signaled by the appearance of a nonzero value of 
IMfcJ for k z — 5?t/12. At a lower temperature T ~ 1.77, a second transition occurs to a lamellar 
phase characterized by a nonzero |MfcJ for k z = n/2 (all other Mk's are zero). 



Our Monte-Carlo results provide strong evidence that the transition from the paramagnetic 
phase is a fluctuation-induced first-order one, in contrast to the Ising ferromagnet with dipolar 
interactions for which the transition is second-order The only (second-order) critical point 
in the temperature- frustration phase diagram is at 7j. Since the transition line approaches 
T c ° from below in a nonanalytic fashion as Q K with n ~ 0.25 — 0.35, avoided critical behavior 
can still dominate the physics of weakly frustrated system above the transition line, thereby 
leading to frustration-limited domain structures in this region. Additional simulation studies 
of the paramagnetic state will however be necessary to confirm this prediction. We have also 
shown that the model gives rise to a variety of modulated phases whose characteristics vary 
both with temperature and frustration. The Laboratoire de Physique Theorique des Liquides 
is UMR No 7600 au CNRS. We are grateful to D. Kivelson, S. Kivelson, Z. Nussinov, A. Lesne, 
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Fig. 4. - Spin configurations corresponding to the different phases observed in the Monte-Carlo 
simulations at Q — 0.144 (box size 12 x 12 x 24): (a) disordered (paramagnetic) phase, (b) modulated 
phase with averaged half-period of 2.4 lattice units; (c) lamellar structure with a half-period of 2 (this 
is also the ground state). Up and down spins are represented by diamonds and dots, respectively. 



and J.-M. Victor for many stimulating discussions 
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